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Abstract 

We derive an integral expression for the plane- wave expansion of the time- varying (nonstation- 
ary) random field inside a mode-stirred reverberation chamber. It is shown that this expansion is 
a so-called oscillatory process, whose kernel can be expressed explicitly in closed form. The effect 
of nonstationarity is a modulation of the spectral field on a time scale that is a function of the 
cavity relaxation time. It is also shown how the contribution by a nonzero initial value of the field 
can be incorporated into the expansion. The results are extended to a special class of second-order 
processes, relevant to the perception of a mode-stirred reverberation field by a device under test 
with a first-order (relaxation-type) frequency response. 



1 Introduction 

In recent years, complex scenarios and applications in electromagnetism have sparked a renewed in- 
terest in the theory and characterization of random electromagnetic (EM) fields in static and dynamic 
environments. Examples include characterization of the impedance and radiation efficiency of electri- 
cally small antennas P], multi-terminal wireless communications systems in multi-path high- mobility 
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propagation environments [2], dynamic atmospheric propagation effects [3] on the group delay for 
satelhte-based navigation systems, radio-frequency and optical speckle and scintillation associated 
with interaction of EM waves with rough surfaces or random media [1], parametrically varying mi- 
crowave billiards and cavities [5]-[7], electromagnetic compatibility (EMC) including testing of emis- 
sions or immunity of electronic equipment to irradiating EM fields [8], etc. Continuously evolving 
multi-scattering electromagnetic environments (EMEs) are ubiquitous but difficult to characterize 
and analyze accurately, even using full-wave numerical simulation methods, because of their large size 
relative to the wavelength and consequent extreme sensivity to small perturbations or configurational 
uncertainties, and because of complicated transient effects occurring in the propagating signals that 
result from dynamically changing boundary conditions, especially in resonant multi-scatter environ- 
ments [7]-[8]. Because of the spatial and/or temporal fluctuations of the wave properties in such an 
EME, their characteristics can often be considered as quasi-random and can be efficiently modelled 
using stochastic methods. From a physics point of view, the task is then to characterize the statistical 
parameters and distributions of the random field in terms of deterministic physical quantities of the 
problem of interest. 

A canonical EME for generating statistically random vector fields is the mode-tuned or mode- 
stirred reverberation chamber (MT/MSRC) [9], [10]. It consists of an electrically large cavity furnished 
with a mechanical stirring mechanism to alter the boundary and/or excitation conditions (e.g., a large 
reflective paddle wheel or moving walls) or a varying (hopping or swept) excitation frequency or noise 
source for mixing the modal fields. Such cavities are important devices for studing the temporal, 
spatial, structural, and polarization properties of local instantaneous scalar and vector wave fields, as 
well as characteristics of their spatio-temporal coherence (correlation) properties. 

For a static or quasi-static MT/MSRC, the assumption of wide-sense quasi-stationarity of the 
field is an excellent and often-made approximation. However, when the rate of change of the EME 
(expressed by the correlation time in space or time) is relatively fast compared to the characteristic 
time or length scale of the field in the corresponding static EME (e.g., frequency, relaxation time), 
this assumption breaks down [11]. In this case, more sophisticated field models and methods of 
analysis are necessary. The purpose of the present paper is to develop a framework for the theoretical 
characterization of such nonstationary fields. 

The issue bears resemblance to the filtered observation of dynamic scenes through natural vision. 
The limited image resolution time of the human naked eye (~ 30 ms) and processing time of vision 
by the brain (~ 150 ms) limits the amount and capacity of visual information that can be faithfully 
captured in real time. In particular, ultra-short-term transient natural phenomena may therefore be 
perceived as distorted (blurred) or even completely masked. This is a manifestation of (nonlinear) 
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nonstationary filtering of the actual scene. Using high-speed cameras, the image frames and dynam- 
ics of the hidden reality can be reconstructed. Parenthetically, compound eyes allow for detecting 
transients more rapidly, through high temporal resolution of flicker across the ommatidia, but at the 
expense of poorer spatial resolution compared to single-chambered eyes [12j . |13j . 

2 Methods for expansion of vector EM fields in dynamic complex 
confined environments 

For the representation of fields inside a MSRC or in a complex EME exhibiting multiple scattering 
or multipath fading, two methods have found application. Under idealized conditions, both methods 
should lead to equivalent results, but each has its particular advantages and disadvantages. Modal ex- 
pansions (MEs) consider the EME as a multi-mode system characterized by a set of vector eigenmodes 
'^mnp discrete resonance frequencies (eigenfrequencies) ujmnp at time t. These eigenmodes form a 
set of basis functions for expanding the local stationary interior field, with expansion coefficients ttmnp* 
In a dynamic EME, the eigenmodes and expansion coefficients both depend on time ^14j : 



MEs are exact, in the sense of inherently incorporating the EM boundary conditions and cavity 
geometry into the description. Since modes are nonlocal, the description inherently includes the 
spatial inhomogeneity of the field caused by the vicinity of boundaries and objects. Since a ME does 
not presume quasi-randomness of the field, it can be used in principle at any frequency, although 
the increasing modal count and modal overlap (modal coupling) in high-Q chambers for increasing 
frequencies may make their implementation prohibitive at very high frequencies. A greater burden 
is the fact that realistic chambers often exhibit complex-shaped time-varying boundaries (diffractors, 
e.g., corrugations and mode stirrers), which hampers the practical implementation of MEs that require 
the calculation and/or functional form of the eigenmodes. Furthermore, a correct description requires 
the cavity to be described as a partially opened system P^, due to the transmitter and receiver 
introducing apertures to the unbounded exterior region of the cavity. This requires mixed (Robin) 
boundary conditions and complicates the analysis and eigenmode characterization. 

In a second method, angular spectral plane-wave expansions (ASPWEs) have been used for deter- 
ministic |16j . [17j and random [18], [19], |20) fields. In ASPWEs, the local field with a suppressed 
exp(ja;t) time dependence for its carrier is expanded as a discrete sum (for bounded closed finite-sized 
cavities) or integral (for unbounded opened regions supporting quasi-random fields) of an angular 
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spectrum of plane waves [16], [21], [22], [23] : 



where = 27r sr or 47r sr is the sohd angle associated with a half-space or entire space for the 
incident field, respectively. Unlike MEs, ASPWEs provide characterizations of the locaU field, thus 
avoiding the need for determining eigenmodes explicitly. They lend themselves easily to statistical 
characterization of this local field and, when extented to spatial correlation functions, to nonlocality 
within a finite region. With the aid of proposed axioms for the average and covariance of the amplitude 
spectrum [3], [21], [25] and spectral density [2lj of the wave components, explicit calculation of average 
and covariance of the expanded random field is possible. By a posteriori imposing EM boundary 
conditions, the ASP WE even enables a complete derivation of probability density functions based on 
the statistical moments |23] . [26j . 

In both MEs and ASPWEs, integral representations are used as approximations to discrete sums in 
cavities that are sufficiently large relative to the wavelength. Strictly, both methods are applicable only 
for static or quasi-statically perturbed cavities, for which the concept of modes and their decomposition 
into standing plane waves are properly defined. When boundaries are moved, complications arise 
relating to nonstationary transient fields and dynamics of eigenmodes, as well as the very concept and 
definition of eigenmodes in such circumstances. Note that, in order for the effect of this motion to 
be significant, this does not necessitate relativisitic speeds, unlike in unbounded EMEs: the relevant 
time scales are the ratios of the correlation time of the mode-stirring process to the decay constant 
of the cavity (at a given frequency), to the response time of the antenna or to the characteristic time 
(e.g., cycle period) of the device under test (DUT), and of course to the period of the source wave 
itself [27|. For example, it is well known that even small Doppler shifts (cf. Sec. 13. ip . frequency 
changes of the order of a few hundred Hz at 1 GHz may cause degradation of performance of wireless 
communication systems, because they give rise to considerable jitter and distortion in impulse radio 
and other wideband signals. 

In this paper, we develop an extension of the ASPWE for steady-state fields to nonstationary 
fields. The analysis is for temporal nonstationarity of scalar fields; corresponding results for spatial 
nonstationarity and vector fields follow mutatis mutandis. 



^If the field is statistically homogeneous (spatially uniform), then the characterization applies of course throughout 
the volume. 
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3 Reverberation process with stirring slip 



3.1 Spectral expansion 

A model of a slipping random field based on a random- walk model was presented in [28] , to which we 
refer for details. Here we briefly review and summarize the main results of this model. 

We consider the evolution of a nonstationary detected statistical mean field Y{t), resulting from 
a stationary source field X{t) (Fig. [T]). Here, X{t) is taken to be a periodic field with deterministic 
waveform, as emitted inside a static cavity. (By extension, a wide-sense stationary noise could be 
treated.) In the application to EMC in Sec. 13.51 X{t) and Y{t) are modulated complex RF fields. 
The field Y{t) is produced by random step changes 5X{t) induced by a change in boundary conditions 
and weighted in the mean by a relaxation time, i.e., it is governed by the first-order Langevin-Ito 
stochastic differential equation (SDE): 

+ (3) 

at T T 

This equation of motion for the random field spans a variety of cases of field dynamics across a time 
interval At, ranging from instantaneous response {At/r +oo) to Brownian motion (At/r 0). 
Within each infinitesimal time increment 5t, Y(t) evolves from its initial value at to in direction 
of its asymptotic steady-state value that would be reached for At — )• +oo. In Fig. [H Y^^\ti) for 
ti G [toi^i] represents the mean outgoing (fading) field for tj, i.e., represented by the set of plane- 
wave components associated with the "previous" state of the system and evaluated at t = tj_. This 
"previous" field decays at a rate governed by the decay time r for each plane- wave component. By 
contrast, Y^^\ti) at any ti G [toi^i] represents the mean incoming (emerging) field at ti, i.e., associated 
with the set of plane- wave components for the "next" state, as evaluated at t = The fact that the 
outgoing field requires a nonzero time to vanish gives rise to a "slipping" field, i.e., a process with a 
gradually fading memory, such that Y{t) is a nonuniformly weighted mixture of both contributions, 
viz., y(to+At) = y(i)(to+At)+c(At)y(°)(to) withO < |c(At)| < 1 (cf. [ig for a detailed classification 
and discussion). 

The fact that this evolution can be characterized, in the mean, by a single time constant r governing 
the rate of change of the instantaneous mean envelope of the transient field (mean effective relaxation 
time r) has been amply supported by numerical and practical experiments [9], [29], |30j . In actuality, 
each component of the field, whether angular-spectral or modal, has its own particular relaxation 
time constant associated with it. However, their values are relatively close within a narrow band of 
frequencies, in the sense that 



'mnp 'm'n'p' 
Tmnv 
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« 1 (4) 




Figure 1: (color online) Discretized sample-and-hold input process X{t) and resulting weighted mean field 
Y{t) = Y'-°^t) for a first-order process, with time- independent average relaxation time r = 0.3 St 

(nonlinear weighting). The actual continuous-time stirring process is obtained as the limit 6t — >■ 0. 



for any pair of modes (or standing plane waves with oppositely directed wave vectors) {m,n,p) and 
{m',n',p') located within this frequency band. Furthermore, strong modal coupling, both in time 
and frequency, causes equalization of their different values Tmnp{fmnp-,U)- For these reasons, a one- 
parameter description is often sufficiently adequate, at least in a mean-square sense, i.e., when the 
field arises as some weighted average, as observed in a coupled multi-modal incoherent system. In any 
case, extensions of the model are possible. For example, a second-order model governing the response 
to step transitions would incorporate not only the average rate of Y{t) approaching its steady-state 
value, but also the average level of overshoot of Y{t) above steady-state during such transitions. 

In [281 Eq. (9)], the resultant Y{t) in a discretized interval was expressed iteratively, as 

a function of X{ti). For the present purpose, it is more beneficial to use the following equivalent 
closed-form expression for Y{t) 



Y(tm ^ tfri+l J 



Y{to) exp 

1 — exp 
-|- 1 — exp 



t-to 



T 



+ exp 



St 

T 

t - t 



n m— 1 



i=0 



T 



T 



T 



(5) 



for m > 1, where tm = to + mSt. For m = 0, the same expression holds, but without the second term. 
Neither Y(t) nor its increments 5Ym{t) = Y{tm+i) — Y^tm) are in general stationary, because both 
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depend on t^. Since we shall further be taking the limit At/r — )• and because we shall be mainly 
interested in late times of observation such that (t — tQ)/5t ^ 1, we restrict further evaluation to the 
discrete time instances t = tm+ii when Eq. ([5]) simplifies to 

N / (m + l)5t\ ( m5t\ 
Y{tm+i) = Y{to)expi-^ _^j+expf — j 

f:exp(+'^)x{U) (6) 



1 — exp 



i=0 



Upon taking the limit 6t/T in Eq. ([6]), the discrete sum reduces to an integral from to to tm, and 
we obtain 

^ ^ vu^ ( tm+i-to \ , , l-exp(-f) 

Y{tm) Y{to)exp + hm 



st^o 5t 

X I exp ( -T^^ \ x{s)ds (7) 



tm, / / Q 



to 



At this point, we introduce the Fourier-Stieltjes spectral expansion of the stationary source field 
X{t) [3T], given by 



X{s) = / exp(ja;s)dZx(u;) (8) 

J —oo 

in which Zx{oj) for an ideal mode-tuned field is assigned the following properties [H], [3T] : 

{dZx{i^)) = 0, (9) 
(dZx(6^i)dZ^(^2)) = -a;2)dFx(6^i,2), (10) 

i.e., Zxi;^) is orthogonal. [Note that Zxioj) represents the accumulated (integrated) spectrum of X, 
not to be confused with its differential, i.e., the Fourier spectrum.] With the substitution u = {tm—s)/T 
and letting tm+i tm = t in Eq. d?]), we arrive after some manipulation at 



/ + _ +„\ r+co 

Y{t) = y(to)exp(^ exp(jL^t) 



X \ \ . " ^ dZxiuj). (11) 

1 +jwr 



Turning attention to the first term in Eq. (jlip . and using general properties of the Fourier transform 
and the Dirac delta distribution [321 Eqs. (2.24) and (2.40)] for X{to) = X{t)6{t - to), this term can 
be assigned a spectral expansion, viz., 

X (to) exp * *° 



-^-^(^^^rexpMt-to)]^^ (12) 
-oo 1+jwr 
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+ 00 



+ 00 



X{uj') exp [iuj'{t - to)] du' 



^re.p[Mt-to)]^^_ (13) 
1 + jwr 

Note that X(to) = ^(^o) because nonstationarity has not yet manifested itself at the start time to, 
whence the output process can then be completely identified with the input process at this instance. 
With an interchange of the order of integration, this yields the general result 

Y{to)exp(^-^—^^ = e{u-t)eM}^t)dZx{i^) (14) 

where 

fl/ , . , , /•+°° exp [^J{t - to)] ^ , 

6'(a;;t) = rexp(-ja;to) / , , . . dw . (15) 

J-oo l+jw'r 

In the particular case where we impose the initial condition i^(to) = 0) for simplicity, only the second 
term in Eq. (Ilip remains and can be represented as an oscillatory process |33j . i.e., 

/ + CJO 
0(t;a;,r)exp(ja;t)dZx(w), (16) 
-OO 

whose kernel 

, A l-exp[-(l+ja;r)^] 
4>{t; u:^ r) = (17) 

is a complex amplitude modulation function expressed in explicit, i.e., closed form. It is seen from Eq. 
()16p that the effect of nonstationarity on the plane-wave expansion is a modulation of the spectral 
field on a time scale that is a function of the cavity relaxation time. 

Fig. [2] shows the time evolution of Eq. (fT7|) at selected values of r = Q/w for an overmoded cavity, 
for narrowband operation at an excitation frequency = 27r x 10^ rad/s. It is seen that larger values 
of r (and, hence, Q) result in longer transition times before (j) reaches a regime, as is intuitively clear. 
The strongest variability of (j){t) occurs when t ~ 27r/a;, as expected. 

As is well known, oscillatory processes provide an evolutionary spectral representation for a special 
class of nonstationary processes [33, Eq. (3.9)], [37]. They apply when nonstationarity is weak (i.e., 
sufficiently slowly evolving), so that the evolution can be considered as a modulation of the original 
process X{t). 

For ujT and t/r +00, we retrieve from Eq. (jlip the classical spectral expansion for 
temporally stationary (but spatially homogeneous as well as inhomogeneous) fields, viz., 

/+00 
exp(ju;t)dZx(^^) =X(t). (18) 
-00 

On the other hand, for r ^ l/w we obtain a sinc(ci;t)-type modulation, in view of the symmetric 
integration limits. 
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3.2 Statistics 



The ensemble mean value of Y(t) is 



{Y{t))=Y{to)exp 



t-to 



(19) 



because (dZx((^)) = 0. The covariance function follows from Eq. (Illh as 



{Yiti)Y*ih)) 
= |y(to)|'exp 



tl + t2 — 2to 



+ 



exp [iLojti - t2)] 
l+a;2r2 

tl + t2 — 2to 



X <1 + exp 



T 



exp - *2)] 



■ exp I - ^^ _ ^° ) exp [-jw (tl - to)] 



exp 



T 

^2 — ^0 

r 



exp [jw (t2 - to)] \ dFx{uj). 



(20) 



This general result is important in the investigation of the effect of nonstationarity on the transfor- 
mation of correlation characteristics and effective number of degrees of freedom. For ti = t2 = t, the 
variance follows as 



4(t) = {\Y{t)\')-\{Y{t))\' 

r . .,-1 r f t-to 

1 — 2 exp 



X COS [uj{t — to)] + exp 



2(t-to) 



dFx(w), 



(21) 



which describes the evolution of fluctuation levels under nonstationarity. We now consider two im- 
portant special cases for the input process. 



3.2.1 Ideal vi^hite noise 

If X{t) is ideal white noise, i.e., dFx{uj) = fxioj)duj = a^duj, then from Eqs. (fT9]) . (pTj) . and [ [38] 
Eq. (3.723.2)], 



VTfT 



X 



1 — exp 



2(t - to) 



(22) 



In the limits t/r — )■ and t/r — > -|-cxo, we retrieve Uy 2TTa\{t — to)/''"^ and ay vra^/r, 
respectively, as expected. 
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3.2.2 Ornstein-Uhlenbeck process 



If X{t) is a first-order process tliat is wide-sense stationary and exponentially correlated with finite 
correlation time T, i.e., 



px{t - to) = exp 



t-tp 
T 



(23) 



so that dFxiuj) = /x(^)da; = 2a^T /[vr (l + ^27-2)]^^^^ ^i^^^ gq. (HH) reduces to 



TT 



1 + exp 



+00 



da; 



(l + a;2r2)(l + a;27^) 



vr 

+00 



exp 



- to) 



cos [a;(t — to)] do; 



(24) 



3 (l + a;2r2)(l + a;27-2)- 

The first integral in Eq. ^ is easily calculated with the aid of [l38j, Eqs. (2.124.1) and (2.161.1)], 
yielding vr/(r + T), i.e., providing a positive contribution irrespective of whether r < T or r > T. 
The second integral in Eq. (j24p is obtained by contour integration of 

exp [-z{t - to)] 



(l-z2r2)(l-z27-2) 

across the left half of the complex z-plane (Fig. [3j). The residues in zi = 1/r and 22 = 1/T are 

rexp(-^: 



(25) 



i2i 



2(r2-r2) 



and R2 



rexp(-^) 



respectively, whence 



aUt) 



MI 
r + r 



1 — exp 



2(T2-r2) ' 
2{t - to) 



(26) 



(27) 



Again, this result holds irrespective of whether r < T or r > T. It is verified that for T/t <^ 1, Eq. 
dJZD reduces to Eq. ([2^ with = 2(TQT/vr. For (t — to)/T <^ 1, we retrieve the linear dependence 
as in the case of ideal white noise under the same condition. Furthermore, for T/r — > 0, it is verified 
that Eq. ^ approaches Eq. 1^ with 2<toT = Traj^. 



3.3 Instantaneous energy density 

The evolution of the nonstationary energy density is of fundamental importance for incoherent fields. 
When considering the instantaneous field as pure functions of time (i.e., not as an analytic function 
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or phasor), the following definition for the associated energy density applies [391 P- 14], 



We{t) 



We{t')dt' 



(28) 



where the latter equality holds for a homogeneous time-invariant medium \D_{t) = e(t) *£_{t) = e£_{t)], 
as we shall further assume. 

The previous analysis for a mean complex field Y{t) holds exactly if Y{t) instead symbolizes the 
mean intensity of a complex field. If we maintain the original definition, however, then for anharmonic 
modulated fields that are being perturbed relatively slowly with respect to a central frequency w, we 
can use the Gabor analytic field representation Y{t) = y{t) exp(— jwt) to write 

d[eY%t)] 



Weit)=Yit) 



dt 



(29) 



Therefore, for dispersionless e, the quantity We{t) for analytic fields is proportional to the field intensity 
|y(t)p. The general expression of We(t) for nonstationary random fields is derived in the Appendix 
and is given by ([^T|) -(|59 p . 

For the magnetic energy density Wh, similar expressions follow by replacing e by /x, with y{t) now 
representing the magnetic analytic field. 



3.4 Structure function 

Random field increments provide a transition between stationary systems and fully-developed nonsta- 
tionary systems (Bachelier-Einstein- Wiener-Levy processes). They may also be useful for describing 
undermoded fields, in which small rapid variations "ride on top" of a slowly varying mean value. 
Their second-order properties are characterized in general by the structure function, which describes 
the covariance of field increments [3l[25]. It follows upon substitution of Eq. ([6]) as 

D{m, n) 

^ ([y(to + rriM) - y(to)] [Y{to + nAt) - y(to)]*) 



\yito)? 



1 — exp 



1 — exp 



tn — to 



exp (iujtm-i) - exp 
exp (jwin-i) - exp 



-oo 
T 

tn-1 
T 



1 + {ujtY 



dFxiuj). 



(30) 
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In writing the second term in Eq. (j30p . we made use of the fact that tm-i — im when At /to — )• 0, and 
Ukewise for tn_i. For m = n, Eq. ([30]) reduces to 



which is formahy equivalent to ()2ip . For example, for an Ornstein-Uhlenbeck process, ()3ip becomes 



3.5 Example: EMC immunity testing using nonstationary mode- stirred vs. quasi- 
stationary mode-tuned reverberation fields 

The above analysis can be applied to determine the response of a linear DUT that is characterized by 
an impulse response function h{t) or a frequency characteristic H{uj) and which has been placed inside 
a mode-stirred reverberation chamber. At any location inside the cavity, the local field consists of 
multipath reflections, i.e., plane waves arriving from isotropically distributed directions with randomly 
polarized direction of polarization and uniformly distributed phases. Within a time interval 5t during 
a mode stirring process, each plane wave undergoes a transition of its directions of arrival, direction of 
polarization and absolute phase. The transition of each parameter can be represented by a trajectory 
in state space. The actual (i.e., physical) rate of fluctuation during this transition is a function of the 
rate of change of cavity perturbation and degree of chaoticity of the cavity. The perceived rate is, in 
addition, a function of the characteristic time (response time) of the DUT to changes of the excitation 



In general, the specific form of H{oj) for an DUT is often unknown. Here, as in [TT|, we simply 
consider a class of DUTs and investigate to what extent a mode-stirred (i.e., continuously varying) vs. 
mode-tuned (i.e., discrete, stepwise changing) excitation field has a difference in effect on the response 
of the DUT. Like the excitation, this response is a random function, hence comparison between both 
cases requires determination of one or more statistical metrics, such as the mean, r/%-confidence 
intervals, statistics of the maximum-to-mean ratio, upward threshold crossing frequency, excursion 
length, etc. 

We shall consider the case of an Ornstein-Uhlenbeck process, as a canonical case, for which explicit 
results can be calculated and from which results for white noise X{t) will follow as a special case. Since 




(31) 




(32) 



field. 



12 



the (real) power spectrum is 



dFxiuj) = fxiuj)duj 



Fxfi 



1 + ^^27-2 

where Fxfl = ZxflZJ^Q, the (complex) amplitude spectrum is 



duj, 



Hence, 



dZx{uj) = zx{uj)duj = Y^^^Tj^duj. 



Y{t) = y(to)exp(-^^) +£^exp(jwt) 



(33) 



(34) 



l-exp[-(l+jL^T) 



Zxflduj. 



(35) 



(l+ja;r)(l+jwr) 

Then the "output" field W{t) of the DUT due to a mode-stirred excitation field, i.e., as measured or 
perceived at one of its test ports, a critical internal component, etc., is then |33l Sec. 6] 



W{t) 



-00 
+00 

-00 
+00 



h{u)Y{t — u) exp [— ja;o(t — u)] du 
Huj+uo (w; t)A{uj + Wo; t) exp(jtjt)dZx(a; + loq) 



(36) 
(37) 



H^+u)o (^;*) exp(ja;t) 
y(to) exp ( - ^ — ^ ) 5{uj + Wo) 



l-exp{-[l+j(w + wo)r]*^} ^ 

+ r. . ■/ : N nr. . ■/ : TTTT'^X.O 



[l+j(w + wo)r] [l+j(w + a;o)7l 
where ojq is any constant frequency and 

Hx{uj;t) -- 

at A = w + ujQ. 



duj 



^ A{\;t-u) . 



(38) 



(39) 



4 Second-order processes 

The response of an overmoded system can only be characterized approximately and in the mean, 
as a first-order system. In reality, the large number of participating modes (i.e., those within the 
instantaneous resonance band of the cavity) and their strong intermodal coupling causes the response 
to be more irregular, requiring a higher-order description. To investigate the effect of increasing the 
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order, we here analyze a second-order system. Because of the resonant nature of the modes, this is 
expected to provide more accurate results. Any higher-order process then follows from cascading of 
first- and second-order systems. 

For the second-order system characterized by 



^ + 2C-^+.lY = K,.l 



(40) 



the step response is the inverse Laplace transform of Y{t) = Kqoj'^/[s (s^ -|- 2^5 -|- oJq)]. The residues 
of the simple poles sq = and si 2 = — C ^ j \/^o ~ lead to the solution 



Y{t) 



Kail 



X COS 



ujQ exp {-(t) 



C"^ t — sin 



UJQ 



(41) 



for < ( < ujQ [For overdamped regime ((" > ujq), the step response involves hyperbolic functions.] 

In principle, the previous analysis of the first-order model can be repeated. To this end, Eq. ([5]) 
could be written as 



Y{t.in < t < tm+l) 

= Y{to)(P"'{At)cl){t-t 



m—l 



\- Ct>{t - tm) 

[l-0(At)]5:r-i-(At)X(t,) 

i=0 

[l-(t>{t-tm)]X{tm) 



(42) 



now with 



ujQ exp {-Ct ) 



cos 



i-Ct- sin-1 



(43) 



in view of a;„ ~ w, t„ ~ r. To obtain the spectral expansion, it is advantageous to consider the Gabor 
analytic signal representation Y{t) = Y'{t) + iY"{t). The step response is then 



Y{t) 



X exp 



ujQ exp (-Ct) 



iiJcoi-et-sm-' 



Wo 



(44) 



in Eq. (jHJ, this approach does not lend itself to 



Because of the presence of the factor 1/ 
express Y{t) easily as an oscillatory or related process, in the spirit of Eqs. (|16|) -(jl7 p for the first-order 
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system. Therefore, we write Eq. (j40p instead as a system of two coupled first-order equations: 

+ iy(t) = ^x{t) 



dy(£) 

dt 



n ' 



(45) 



dZ{t) 
dt 



thereby assuming that 



Ti < r2 < 1 



(46) 



where the response is now given by Z{u}), with Y{uj) denoting an auxihary intermediate process. For 
an overdamped second-order system {( > wg), such a decomposition is always possible. The process 
Y{t) has a physical meaning: for a first-order system responding to a nonstationary field which, in 
turn, is governed by a first-order process, the excitation of this system is given by Y{t). The spectral 
representation follows by replacing X{s) in Eq. ([7]), as given by Eq. ([8]) and Y{s) given by Eq. ([TT]) : 



Z{to) exp 



t-tp 

T2 



^ Y{to) 





_to\ 


. \r2 


TlJ. 



X [exp (-2Cto) - exp (-2(1)] 



+ 00 



exp (jwt) - exp (-^) exp (jwto 



(1 +ja;ri) (1 +ja;r2) 



^ exp (ja;^o) n 
1 + jwTi r2 
xdZx(a;). 



exp 



-2Ct + 



exp 



T2 



(47) 



For the case where to = 0, Y{tQ) 

Z{t) = 



Z{to) = 0, this expression reduces to 
hoo Jexp(ja;t) -exp(-:;|) 1 



n 



(1 + jwri) (1 + ju;r2) 
t 



exp {—2(t) — exp 



T2 



1 + jwri r2 
dZx(a;). 



(48) 



For an Ornstein-Uhlenbeck process, contour integration yields 



Z{t) 



-J 



2aiT' 



T 



T + T [T +T2 

n 



exp ( -— ) + exp ( - — 



T2 



exp (-2Ct) 



t 

f 

exp [ - — 
r2 



(49) 



Note that, even though the original ojq and ( have been replaced by ri and r2, exhibiting a one-to-one 
relationship, the process of replacing the second-order SDE by two coupled first-order SDEs requires 
for one time constant to be much smaller than the other one, in order to maintain the validity of the 
Langevin-Ito formulation. 
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5 Conclusions 



Angular spectral plane-wave expansions provide a framework for the stochastic characterization of 
dynamic random EM fields that result from moving boundaries, spatial scanning, and/or changing 
excitation frequency with consequent changes in the pattern of participating eigenmodes. It was shown 
that the spectral expansion can be expressed as an oscillatory process Eq. (jl6p with kernel Eq. (jl7p 
(or, more generally, as Eqs. pT|) . (fT5|) ). The fact that the expansion can still be represented by 
an oscillatory process (which is usually associated with narrowband, i.e., slowly modulating (quasi- 
stationary) processes, but now with no such limitation on the rate of fluctuations) can be attributed to 
the linearity of the system. The results for highly overmoded resonances apply in the mean (incoherent 
superposition), i.e., for the collection of cavity modes participating in forming the instantaneous field 
at any one time. 

A useful extension would be the "microscopic" plane- wave spectral expansion for a single second- 
order resonant system (cavity mode), without any limiting condition on the separation of the two 
time constants involved in decomposing into two first-order processes. Such a complete second-order 
characterization could form the basis for a coherent superposition of excited cavity modes. 
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Figure 2: (color online) Modulating function (j>{t;uj,T) of the kernel in Eq. ([T6|) for the spectral plane-wave 
representation of nonstationary mode-stirred fields, for selected values of t at uj /{2tt) — 1 GHz: (a) magnitude; 
(b) phase. 
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Im(z) 




Re(z) 



Figure 3: Contour of integration for Eq. ((25|) in the complex z-plane {t > T). 
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Appendix: Instantaneous energy density 



Here, we derive the general expression for the instantaneous energy of a nonstationary field. Substi- 
tution of the general expression Eq. ([TT]) for Y{t) into Eq. ([29|) yields 



We{t') 



T 
T 



2{t' - to) 



■ exp 



t'-to\f+'^ l-exp[-(l+ja;r)^^ 
exp [ / exp (jwt ) — 



T 



r 



, eY{to) ( t'-to\ , . {-j^T + exp 

-\ exp / exp (—]UJt * 



dZx{uj] 



g r+oc r+oc 1 — exp — (1 + JcJit) — 

+- / / exp [j {ui - UJ2) t'] -TV- ~ 



1 -iuT 

t'-to' 



-jw2T + exp 



:i-j-2r)^^]} 



1 - ja;2r 

Integration of We{t') yields the instantaneous We{t) as 



dZx{uJi)dZ*x{L02). 



We{t) = f wS')<^t' = We, (t) + We, (t) + We, (t) + We, (t) 
Jto 



e\Yito) 



with 

WeAt) = 

We,{t) = eY*{to: 
eY{to) 



1 — exp 



(50) 



(51) 



2(t-to) 



exp (jcjto) 



We,{t) 
WeAt) 

and 



-00 1 + ^V2 

2{t - to 



exp 



(l-iUT) 



t-to 





r 2{t-to)i 




— exp 


T 


} 



dZxiuj) 



1 — exp 



+00 



exp (jwto) 



1 - jwr 



dZ*x{io) 



\ jaJrexp(-ja;to) 



-00 

hoo i"+oo 



(l+jwr) 



t-to 



dZ*x{u; 



(52) 
(53) 

(54) 



00 J —00 



Il{t;UJi,UJ2) + l2{t;UJi,UJ2) + l3{t;UJi,UJ2) + l4{t;UJi,UJ2) 

(1 +jLJir)(l - jcj2r) 



dZx(^^i)dZ^(w2) (55) 



Ii(t;cji,cj2) 
l2(t; wi,W2) 
l3(t; wi,tJ2) 
l4(t; u;i,tJ2) 



W2 



■ exp [j(a;i - a;2)to] {1 - exp [j(a;i - uj2){t - to)]} 



U)i - UJ2 

- — ] exp ( - — ) exp [j(a;i - ^^2)^0] "i 1 - exp 

j^2T 



T 



(1 - jWlT 



t-to 



1 + ja;2r 



exp [j(u;i -uj2)to]\l- exp 



[l+iuJ2T] 



t-to 



\ \~r' ' ^■^^^^ ~ ^^''^"^ <{ 1 - exp 



r 

2(t-to) 



(56) 
(57) 
(58) 
(59) 
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For oji = U2 = uj, 



hit; to) 



iuj{t - to) 

exp I M -'^ ~ 6xp 



r 

1 — exp 



1 - jwr 

1 +iuJT 

1 / 2to 



(l-jcur) 



t-to 



■(l+ja;r) 
1 — exp 



t-to 



T 

2{t - to) 



Considering that {dZ{u)) = and {dZ [ojifdZ* {uj2)) = 5{oJi — 0J2)dFx{(jj), we have that 



(We,) = {We3)=0 

'■+°° h{t;u) + h{t;uj) + h{t;uj) + h{t;uj) 



{We,) = ef 

J — ( 



1 + OJ^T 



2^2 



dFxiu;). 



(60) 
(61) 

(62) 

(63) 



(64) 
(65) 
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